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Abstract 

We extend a recent chiral approach to nuclear matter by including the most general 
(momentum-independent) NN-contact interaction. Iterating this two-parameter contact- 
vertex with itself and with one-pion exchange the emerging energy per particle exhausts all 
terms possible up-to-and-including fourth order in the small momentum expansion. Two 
(isospin-dependent) cut-offs Aq^i are introduced to regularize the (linear) divergences of 
some three-loop in-medium diagrams. The equation of state of pure neutron matter, 
En{kn), can be reproduced very well up to quite high neutron densities of p„ = 0.5 fm~'^ 
by adjusting the strength of a repulsive nn-contact interaction. Binding and saturation 
of isospin-symmetric nuclear matter is a generic feature of our perturbative calculation. 
Fixing the maximum binding energy per particle to —E{kfQ) = 15.3 MeV we find that any 
possible equilibrium density po lies below p^^^ = 0.191 fm~^. The additional constraint 
from the neutron matter equation of state leads however to a somewhat too low saturation 
density of po = 0.134 fm~^. We also investigate the effects of the NN-contact interaction 
on the complex single-particle potential U{p, kj) + iW{p, kj). We find that the effective 
nucleon mass at the Fermi-surface is bounded from below by M*{kfQ) > 1.4M. This 
property keeps the critical temperature of the liquid-gas phase transition at somewhat too 
high values Tc > 21 MeV. The downward bending of the asymmetry energy A{kf) above 
nuclear matter saturation density is a generic feature of the approximation to fourth order. 
We furthermore investigate the effects of the NN-contact interaction on the (Vp)^-term 
in the nuclear energy density functional £[p,t]. Altogether, there is within this complete 
fourth-order calculation no "magic" set of adjustable short-range parameters with which 
one could reproduce simultaneously and accurately all semi-empirical properties of nuclear 
matter. In particular, the conditions for a good neutron matter equation of state and for 
good single-particle properties are mutually exclusive. 

PACS: 12.38.Bx, 21.65.+f, 24.10.Cn, 31.15.Ew 



1 Introduction 

The present status of the nuclear matter problem is that a quantitatively successful description 
can be achieved, using advanced many-body techniques, in a non-relativistic framework when 
invoking an adjustable three-nucleon force p. Alternative relativistic mean-field approaches, 
including either non-linear terms with adjustable parameters or explicitly density-dependent 
point-couplings, are also widely used for the calculation of nuclear matter properties and of 
finite nuclei |21 El IH E] • At a more basic level, the Dirac-Brueckner method of ref . [H] solves a 
relativistically improved Bethe-Goldstone equation with one-boson exchange NN-interactions. 

^Work supported in part by BMBF, GSI and DFG. 
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In recent years a novel approach to the nuclear matter problem based on effective field 
theory (in particular chiral perturbation theory) has emerged [3 IHl E! • The key element there 
is a separation of long- and short-distance dynamics and an ordering scheme in powers of small 
momenta. At nuclear matter saturation density the Fermi- momentum kfQ and the pion mass 
m^r are comparable scales (/c/o — 2m^), and therefore pions must be included as explicit de- 
grees of freedom in the description of the nuclear many-body dynamics. The contributions 
to the energy per particle E{kf) of isospin-symmetric nuclear matter as they originate from 
chiral pion-nucleon dynamics have been computed up to three-loop order in refs.[7, ^. Both 
calculations include the Ivr-exchange Fock-diagram and the iterated Ivr-exchange Hartree- and 
Fock- diagrams. In ref.|8j irreducible 27r-exchange has also been taken into account and a mo- 
mentum cut-off A has been used to regularize the few divergent parts associated with chiral 
27r-exchange. The resulting cut-off dependent contribution [E{kf) ~ A/cp see eq.(15) in ref.|H]) 
is completely equivalent to that of a zero-range NN-contact interaction. At that point the (ear- 
lier) work of Lutz et al. follows a different strategy. Two zero-range NN-contact interactions 
(acting in the ^Si and ^5*0 NN-states) proportional to the dimensionless parameters Qq + g\/4: 
and Qi + g\/4: are explicitly introduced (see eq.(4) in ref.[Zj). The components proportional 
to g\/4: cancel (order by order) zero-range contributions generated by one-pion exchange. The 
other components proportional to go and gi are understood to subsume all non-perturbative 
short-range NN-dynamics relevant at densities around nuclear matter saturation density po- 

Despite their differences in the treatment of the effective short-range NN-dynamics both 
chiral approaches Ej are able to reproduce correctly several empirical nuclear matter prop- 
erties (saturation density po, binding energy per particle —E{kfo) and compressibility K = 
kjQE"{kfo)) by fine-tuning only one parameter, either the coupling go + gi — 3.23 or the cut-off 
A ~ 0.65 GeV. The fine-tuning of a momentum cut-off A as done in ref.pj introduces some 
degree of non-uniqueness that lies outside chiral perturbation theory. One should, however, 
keep in mind that fine-tuning of some parameter(s) is involved in most calculations of nuclear 
matter. 

Although the conceptual treatment of the short-range NN-interaction by Lutz et al. ^ may 
be very appealing at first sight it suffers (after closer inspection) from severe phenomenologi- 
cal problems. As shown recently in ref.[10j the single-particle properties come out completely 
unrealistic in this scheme. The potential depth of f/(0, fcjo) = — 20MeV is by far too weakly 
attractive while on the other hand the imaginary single-particle potential W{0, kfo) = 51 MeV 
is too large. Most seriously, the total single-particle energy Tkin(p) + U{p,kfQ) does not rise 
monotonically with the nucleon momentum p, thus implying a negative effective nucleon mass 
at the Fermi-surface p = kfo (see Fig. 3 in ref.jTni). In such an abnormal situation there exist 
occupied nucleon- states in the Fermi-sea with total energy higher than the Fermi-energy, indi- 
cating an instability of the system. As a matter of fact the overly strong momentum dependence 
of the nuclear mean-field U{p, kf) in the scheme of Lutz et al. [Zj (leading to dU /dp < —p/M) 
originates from those diagrams in which the (strong and) attractive NN-contact interaction 
(proportional to g^ + gi + g\/'2.) is further iterated. The NN-contact interaction necessarily 
has to be very strong in this scheme since it effectively also has to account for the attraction 
generated dynamically by iterated Ivr-exchange in form of a (regularization dependent) linear di- 
vergence d/ 1. In contrast to that much better single-particle properties have been obtained 
in the chiral approach of ref. [0] employing cut-off regularization and not introducing any explicit 
short-range NN-terms. For example, the resulting potential depth of ?7(0, kf^) = —53 MeV jH] 
is in good agreement with the depth of the empirical nuclear shell model potential or the 
optical model potential ^2j. Furthermore, when extended to finite temperatures this approach 
reproduces the first-order liquid-gas phase transition of isospin-symmetric nuclear matter, how- 
ever, with a somewhat too high value of the critical temperature Tc = 25 MeV ^3] • The reason 
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Figure 1: Additional in- medium diagrams generated by the NN-contact interaction. The filled 
square vertex symbolizes this zero-range NN-contact interaction proportional to 70 and 71. The 
indices 1 = 0,1 label the total isospin of the NN-system. 

for that is a still too strong momentum dependence of the real single-particle potential U {p, kjo) 
near the Fermi-surface p = kfo [9^. 

The purpose of the present paper is to investigate in detail the role of short-range NN-terms 
in addition to pion-exchange dynamics for a variety of nuclear matter observables. Our primary 
interest is to find out whether such adjustable NN-contact terms (of moderate strength) can 
cure certain shortcomings of previous calculations jHl IH which have included only the pion- 
exchange contributions to three-loop order. From the point of view of the expansion in powers 
of small momenta {kf and m^) the combination of pion-exchange and zero-range NN-contact 
interaction exhausts actually all terms possible up-to-and-including fourth order. Alternatively, 
the interacting part of the energy per particle can be understood (within this approximation 
to order 0{kj)) as to result from the Hartree and Fock contributions of a nucleon-nucleon 
T-matrix of the form: 



NN 



^ ^ 7o + 371 7i - 7o ^ 
T"i • r2 H \ Ti 



1"2 



(1) 



which is evaluated in first and second order perturbation theory (i.e. it is also iterated once 
with itself). Note that additional spin-spin terms proportional to ai ■ 02 are redundant in the 
T-matrix eq.(l) as a consequence of the Pauli exclusion principle. The two coefficients, 70 
and 7i, parameterize the strength of the contact interaction in NN-states with total isospin 
/ = 0, 1. In order to have immediately a measure of the (relative) magnitudes involved we 
note that the repulsive contact piece which survives from one-pion exchange in the chiral limit 

= can be mapped onto 70"^^^^ = 7(^'^'' = —1/2. As usual, j-,, = 92.4 MeV denotes the pion 
decay constant and = 135 MeV stands for the (neutral) pion mass. The value qa = 1-3 of 
the nucleon axial-vector coupling constant is linked via the Goldberger-Treiman relation to a 
vrA^A^-coupling constant of q^^n = Qa^/ fir = 13.2 The iteration of the T-matrix eq.(l) to 
second order generates besides finite terms also the (primitive) linearly divergent loop-integral 

dl 1. As done in ref.^j it will be regularized by a momentum cut-off Ao,i which we allow in 
this work to depend on the total isospin J = 0, 1 of the two-nucleon system.^ Such a doubling 
procedure is physically meaningful since cut-off dependent contributions are fully equivalent to 
contributions from (momentum independent) NN-contact vertices. On the practical side four 
free parameters, 70,1 and Ao,i give us more flexibility in adjusting nuclear matter properties 
(within the present perturbative scheme). Of course, the two cut-offs Ao,i should not differ too 
much from each other and they should also take on physically reasonable values. 

^The weighting factors of Aq and Ai are readily found by decomposing the isospin factor of a particular 
diagram into the projection operators on total isospin / = 0, 1: Vi:, = {1 — ti ■ T2)/4 and Vi = {i + fi ■ T2)/4. 
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2 Neutron matter equation of state 



We start the discussion with the equation of state of pure neutron matter. For this system 
the total isospin of two nucleons is restricted to / = 1 and therefore only the strength param- 
eter 7i (see eq.(l)) and the isospin-1 cut-off Ai can come into play. We first write down the 
contributions to the energy per neutron En{kn) generated by the zero-range nn-contact inter- 
action proportional to 71. The left diagram in Fig. 1 obviously leads to a contribution linear 
in the neutron density p„ = /c^/Stt^. Of the same structure are all the cut-off dependent terms 
from iterated diagrams with two medium insertions^. We combine them by introducing the 
dimensionless parameter Fi: 



En{kn) 



Ti = 7i + 7T-7T^(27i - 1 • 



(2) 



The second diagram in Fig. 1 with two medium insertions on parallel nucleon lines gives rise to 
a finite (i.e. cut-off independent) contribution of the form: 



5(87r)3/4 



22u 4(5 + Au^) arctan 2u 

u 



+ - ) ln(l + 4u2) 



Au^ u ) 



(3) 



with u = kn/mTT the ratio of the neutron Fermi-momentum fc„ and the pion mass 171^^. These 
two quantities are treated as small momentum/mass scales in our perturbative calculation. 
The second and third diagram in Fig. 1 with three medium insertions represent Pauli-blo eking 
effects. The corresponding contribution to the energy per neutron reads: 



En{kn) 



(47rA)%3 



7iM 
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l + {u + x) 
1 + {u — x) 

U + X 



In 



U — X 



, (4) 



with the auxiliary functions H = hi{u + xy) — ln{u — xy) and s = xy + a/m^ — + x'^y'^. The 
quantity s has the geometrical meaning of the distance between a point on a sphere of radius 
u and an interior point displaced at a distance x from the center of the sphere. In the same 
geometrical picture y denotes a directional cosine. Note that iterated three-loop diagrams carry 
an energy denominator equal to the difference of nucleon kinetic energies and therefore they 
are proportional to the (large) nucleon mass M = 939 MeV. The additional contributions to 
En{kn) from the (relativistically improved) kinetic energy, from Ivr-exchange and from iterated 
iTT-exchange have been written down in eqs. (32-37) of ref.|H]. In the case of the Ivr-exchange 
contribution (see eq.(33) in ref.jH]) we neglect of course the small relativistic l/M^-correction of 
order 0{kj). Note that the parameter Fi defined in eq.(2) includes also the cut-off dependent 
contribution from iterated Ivr-exchange (see eq.(39) in ref.jH]) through its 71-independent piece. 

The full line in Fig. 2 shows the energy per particle En{kn) of pure neutron matter as a 
function of the neutron density pn = k^/Sn^. The dashed line stems from the many-body 
calculation of the Urbana group ^J. This line should be considered as a representative of the 

■^This is a technical notation for the difference between the in-medium and vacuum nucleon propagator. For 
further details, see sect. 2 in ref.(S|. 
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Figure 2: The energy per particle En{kn) of pure neutron matter versus the neutron density 
pn = ^n/3vr^. The full line is obtained by adjusting the strength parameter of the rm-contact 
interaction to 71 = —0.75 together with a cut-off Ai = 0.48 GeV. The dashed line stems from 
the many-body calculation of the Urbana group . The dotted line gives one half of the kinetic 
energy Ekin{kn)l2 = 3k^j20M. 

host of realistic neutron matter calculations ^SJ ^2] which scatter around it. The full line 
in Fig. 2 has been obtained by adjusting the parameter 71 to the value 71 = —0.75. This 
number translates into a repulsive zero-range nra-contact interaction which is comparable in 
strength to the Ivr-exchange. (Remember that the contact piece of Ivr-exchange maps onto 
^(Itt) _ _]^^2). The other parameter Fi regulating an attractive contribution to En{kn) linear 
in the neutron density pn takes the value Fi = 2.76. The cut-off scale Ai ^ 0.48 GeV which 
parameterizes (via the relation in eq.(2)) this amount of attraction linear in density lies well 
within the physically reasonable range. One observes that the result of realistic neutron matter 
calculation of ref.^] can be reproduced up to quite high neutron densities pn < 0.5 fm~^. In 
particular, the downward bending of En{kn) above pn > 0.2 fm~'^ which showed up in Fig. 8 of 
ref.^ (including only pion-exchanges) is now eliminated. Note that in the scheme of Lutz et 
al. [Zj this downward bending also persisted (see Fig. 6 in ref. [THj). Altogether, it seems that 
most of the dynamics underlying the neutron matter equation of state (at moderate densities 
Pn < 0.5 fm~^) can be covered by pion-exchange and a comparably strong short-distance nn- 
repulsion treated up to three-loop order. At low densities, p„ < 0.2 fm~^, the equation of state 
of pure neutron matter is surprisingly well approximated by just one half of the kinetic energy, 
Ekin{kn)/2 = 3/c^/20M, (see the dotted line in Fig. 2) as emphasized recently in ref . fTT] . 

3 Equation of state of isospin-symmetric nuclear matter 

In this section we turn to the equation of state of isospin-symmetric nuclear matter. In this 
system with proton-to-neutron ratio equal to one both total NN-isospins / = 0, 1 are present 
and therefore all four parameters 70,1 and Ao,i come into play. Following the scheme in sect. 2 we 
first enumerate the contributions to the energy per particle E{kf) generated by the zero-range 
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NN-contact vertex. The complete term linear in the nucleon density p = 2kj-/3iT'^ can again 
be compactly written by introducing a new parameter Fq which subsumes cut-off dependent 
pieces in the 1 = channel: 

The second pion-exchange diagram in Fig. 1 with two medium insertions on parallel nucleon 
lines leads to the following finite (i.e. cut-off independent) contribution: 



> 3(70 + I'MMmt 



22m - - 4(5 + 4^2) arctan2M + ( 7^ + - ) ln(l + 4^^^ 



V4m3 U J 



(6) 



where the meaning of u has changed in this and all following sections to m = kf/m.,^. Diagrams 
with three medium insertions represent Pauli blocking effects in isospin-symmetric nuclear 
matter. Their contribution to the energy per particle reads: 



X (7o + 7i - 70 - 7i)s' + (70 + 71) ln(l + s^)] , (7) 

where the linear/quadratic term in 70,1 obviously belongs the second/third diagram in Fig. 1. 
For further reduction of the occurring double- integral, see eq.(4). The expansion of the energy 
per particle up to order 0{kj) is completed by adding to the terms eqs. (5,6,7) the contri- 
bution from the (relativistically improved) kinetic energy, (static) Ivr-exchange and iterated 
Ivr-exchange written down in eqs. (5-11) of ref.jH]. 

Let us first look at generic properties of the nuclear matter equation of state in our calcu- 
lation. Binding and saturation occurs in a wide range of the three relevant combinations of 
parameters: Fq + Fi, 70 + 71 and 7o + 7i • We fix the negative minimum of the saturation curve 
E{kf) to the value —15.3 MeV [TH] by adjusting at given (70, 71) the strength Fq + Fi of the term 
linear in the density p appropriately. The horizontal position of each minimum determines then 
the corresponding saturation density po(7o,7i) = Po(7i,7o)- It is obviously symmetric under 
the exchange of the two parameters 70 ^ 71. In Fig. 3, we show contours of constant satura- 
tion density po in the (70, 71) -plane. To a good approximation these contours are concentric 
circles around a point in parameter space where the saturation density becomes maximal. We 
find numerically the following approximate parameter dependence of the saturation density (in 
units of fm~'^): 

Po[fm-^] = 0.191 - 0.0388((7o - 0.387)^ + (71 - 0.387)^) . (8) 
Note that the maximal saturation density p^^^ = 0.191 fm~'^ possible in the present perturbative 



calculation lies just at the upper end of the range po = (0.166 ± 0.027) fm"^ quoted in ref.jll 
Interestingly, the maximal saturation density p™^^ = 0.191 fm~^ is reached in the near vicinity 
of the special point 70 = 7i = 1/2. The latter describes the situation in which the (explicitly 
introduced) short-range NN-interaction exactly cancels the contact pieces of the Ivr-exchange. 
The additional constraint 71 = —0.75 taken over from the fit to the neutron matter equation 
of state (see sect. 2) leads (together with 70 = 0.39) however to a somewhat too low saturation 
density of po = 0.134 fm~^. The other parameter Fq = 5.63 needed to keep the depth of the 
saturation curve at —15.3 MeV translates into an isospin-0 cut-off of Aq = 0.55 GeV. Its value 
does not differ much from the isospin-1 cut-off Ai ~ 0.48 GeV needed for the neutron matter 
equation of state and it also lies within the physically acceptable range. Note that momentum 
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Figure 3: Contours of constant nuclear matter saturation density po (in units of fm ^) in 

the (70, 7i)-plane of short-range parameters. The maximal saturation density is p^^^ = 

0.191 fm"^ ~ l.lSpo'^P. The binding energy per particle is fixed to the value — i?(fc/o) = 
15.3 MeV. 



cut-offs around 0.5 GeV have also been used in chiral perturbation theory calculations of NN- 
phase shifts in refs.|2()^ 1^ . In these works the cut-off serves a somewhat different purpose 
of regulating the (non-perturbative) Lippmann-Schwinger equation. However, no comparable 
fine-tuning of A is needed in these works since the cut-off dependence of NN-observables can 
to a large extent be balanced by the strengths of NN-contact terms. 

One may also invert the fitting procedure and extract this way the parameter dependence 
of the maximal binding energy —E{kfQ) per nucleon at fixed saturation density po- For this 
purpose we choose the value po = 0.158 fm~^ as obtained by extrapolation from inelastic electron 
scattering off heavy nuclei [22j. Fig. 4 shows the resulting contours of constant saturation energy 
E{kfQ) in the (70, 7i)-plane of short-range parameters. Again, to a good approximation these 
contours are concentric circles around the point 70 = 7i = 0.38 in parameter space where the 
saturation energy becomes maximal, £'(/c/o)™^^ = —10.7 MeV. 

The full curve in Fig. 5 shows the equation of state of isospin-symmetric nuclear matter cor- 
responding to the parameter set (70 = 0.39, 71 = —0.75, Fq = 5.63, Fi = 2.76) constrained by 
the equation of state of pure neutron matter. The nuclear matter compressibility (related to the 
curvature of the saturation curve at its minimum) comes out as K = kj^E^kfo) = 243 MeV, 
in good agreement with the presently accepted empirical value K = (250 ± 25) MeV fOl 
For comparison, the value K = 253 MeV ^0] of the nuclear matter compressibility is obtained 
if the (explicit) short-range NN-terms are switched off (70 = 71 = 0, Aq = Ai = 0.61 GeV). 
The dashed line in Fig. 5 shows the equation of state E{kf) with maximal saturation density 
pmax _ Q ]^g]^ fm~'^ which is realized by the parameter choice: 7o = 7i = 0.39 and Fq + Fi = 6.46 
(or equivalently Aq = Ai = 0.60 GeV). In this case the nuclear matter compressibility changes 
only slightly to i^' = 250 MeV. The dotted line in Fig. 5 stems from the sophisticated many-body 
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Figure 4: Contours of constant saturation energy Eikfo) (in units of MeV) in the (70, 7i)-plane 
of short-range parameters. The maximal saturation energy is Eikfo)^^^ = —10.7 MeV. The 
nuclear matter saturation density is fixed to the value po = 0.158 fm~^ j22] . 

calculation of the Urbana group jl] using the Argonne V18 NN-potential plus boost corrections 
and three-nucleon interactions. By giving up the constraints from the neutron matter equa- 
tion of state one could also reproduce this curve (for moderate densities p < 0.3 fm~^). The 
corresponding best fit parameters would be 70 = 7i = —0.19 and Aq = Ai = 0.59 GeV. Alto- 
gether, it seems that most of the parameter freedom of our calculation is already exhausted by 
adjusting the nuclear matter saturation point. The compressibility is rather stable and takes 
on reasonable values around K ~ 250 MeV. 



4 Real single-particle potential 

Next, we turn to the real part of the single-particle potential U{p, kf) below the Fermi surface 
{p < kf). As outlined in ref.p, the contributions to the (real) nuclear mean-field U{p,kf) 
can be classified as two-body and three-body potentials. The first diagram in Fig. 1 (combined 
with all cut-off dependent terms from iterated diagrams) leads to a momentum independent 
contribution to the two-body potential of the form: 

(ro + r,)g\k} 
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Figure 5: The energy per particle E{kf) of isospin-symmetric nuclear matter versus the nucleon 
density p = 2/cj/37r^. The full curve corresponds to the parameter set constrained by the 
neutron matter equation of state. The dashed curve shows the case of maximal saturation 
density p^^^ = 0.191 fm~^. The dotted line stems from the many-body calculation of the 
Urbana group [Tj. 



which is just twice the energy per particle written in eq.(5). From the second pion-exchange 
diagram in Fig. 1 one derives a (finite) contribution to the two-body potential of the form: 



(70 + li)9AMmi 
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(10) 



with the abbreviation x = p/rriT^. The second and third diagram in Fig. 1 give each rise to three 
different contributions to the three-body potential. Altogether, they read: 
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(11) 



with the auxiliary function a = ^y ^ a/m^ — -|- ^y'^ and the ratio of parameters r = (7q -t- 
7i)/(7o + 7i)- The real single-particle potential f/(p, kj) is completed by adding to the terms 
eqs. (9,10,11) the contributions from l7r-exchange and iterated Ivr-exchange written down in 
eqs.(8-13) of ref.||9j. Again, the higher order relativistic 1/M^-correction to Ivr-exchange (see 
eq.(8) in ref-jH]) is neglected here for reasons of consistency. 
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Figure 6: The lower curve shows the real part of the single-particle potential U{p,kfo) at 
saturation density kfo = 248 MeV for the parameter set constrained by the equation of state 
of pure neutron matter. The upper curve includes the relativistically improved kinetic energy 
rkin(p) =pV2M-p78MI 



The slope of the real single-particle potential U{p, kf) at the Fermi surface p = kf determines 
the effective nucleon mass (the product of "fc-mass" and "£'-mass" divided by the free mass 
in the nomenclature of ref.jSl]) via a relation M*{kf) = M/(l + B - kj/2M^). For the slope 
parameter B we find the following numerical result (choosing kf = 2m^): 

= -0.239 - 1.245((7o - 0.388)^ + (71 - 0.388)^) . (12) 

p=kf=2m,Tr 

The negative coefficient in front of the quadratic polynomial in 70 and 71 is known analytically 
from Galitskii's calculation |2S1I2E] and it reads in our notation: g\M'^k'j{l — 71n2)/(807r^/^). 
One deduces from eq.(12) that the slope parameter B is negative definite. This implies that 
according to our calculation the effective nucleon mass at the Fermi-surface M*{kfQ) is al- 
ways larger than the free nucleon mass M = 939 MeV. While the iteration of Ivr-exchange 
with the NN-contact vertex (second diagram in Fig. 1 being linear in 70 + 71) can of course 
reduce the effective nucleon mass, there are counteracting quadratic effects from the iteration 
of the NN-contact interaction with itself, which set a limit to this. The lower bound on the 
effective nucleon mass in the present calculation is M*(/cjo) ^ 1.4M. It exceeds considerably 
the value M*(fcjo) — 1.15M found in the self-consistent Brueckner calculation of ref.j^Zj (see 
Fig. 6 therein). Interestingly, the parameter set (70 = 71 = 0.39) leading to the minimal ef- 
fective nucleon mass coincides with that of maximal saturation density p™^^ = 0.191 fm~^. 
Both features occur simultaneously in a situation where the (explicitly introduced) short-range 
NN-interaction cancels almost completely the contact piece of the Ivr-exchange. 

The lower curve in Fig. 6 shows the real part of the single-particle potential f/(p, fc/o) at 
saturation density fc/o = 248 MeV for the parameter set constrained by the equation of state 
of pure neutron matter. The potential depth of f/(0, fc/o) = —27.7 MeV is by far too weakly 
attractive, i.e. at least a factor of two too small. Moreover, a very strong decrease of U{0, kfo) 
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Figure 7: The lower curve shows the real part of the single-particle potential U{p,kfo) at 
kfQ = 279 MeV for the case of maximal saturation density p^^^ = 0.191 fm~^. The upper curve 
includes the relativistically improved kinetic energy T\^in{p) = /2M — p'^/8M^. 

sets in above p = kfo/2 = 124 MeV. The consequence of this behavior is that the total single- 
particle energy Tkin(p)+f/ {p, kfo), shown by the dashed line in Fig. 6, does not rise monotonically 
with the nucleon momentum p. From that perspective the single-particle potential U{p,kfQ) 
presented in Fig. 6 is as unrealistic as the one in the scheme of Lutz et al. [71 ^2] • The 
requirements of a good neutron matter equation of state and good single-particle properties in 
isospin-symmetric nuclear matter cannot be fulfilled simultaneously in a calculation to fourth 
order. A possible remedy of this could be a p^-dependent NN-contact interaction (contributing 
at fifth order in small momenta). 

The lower curve in Fig. 7 shows the real part of the single-particle potential U{p,kfo) for 
the case of maximal saturation density p^^^ = 0.191 fm~'^. The potential depth of U{p, kjo) = 
—65.5 MeV is now close to the value —64 MeV obtained in the calculation of ref. [2Z] using the 
Paris NN-potential. The momentum dependence is also reasonable since its leads to a monotonic 
rise of the total single-particle energy Tkin(p) + ^/o), as it should be. Note that both dashed 
curves in Figs. 6,7 end at the Fermi surface p = kfo with the value E{kfQ) = —15.3 MeV as 
required by the Hugenholtz- Van-Hove theorem [2H]- This serves an important check on our 
(analytical and numerical) calculations. 

5 Imaginary single-particle potential 

Next, we come to the imaginary single-particle potential W{p, kf) as it arises from the iteration 
of the NN-contact interaction with Ivr-exchange and with itself. For p < kf the imaginary 
single-particle potential W{p, kf) is equal to the half-width of nucleon hole states in the Fermi 
sea. As outlined in ref.0 the contributions to W{p,kf) can be classified as two-body, three- 
body and four-body terms. From a technical point of view a u-hodj term involves an integral 
over the product of z/ — 1 Fermi spheres of radius kf. The second and third diagram in Fig. 1 
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give rise to a two-body term of the form: 



(r - 1) 2u'x^ - 3uM + 4 + Uu^ - 

5 / 3 



+ — (3x^ — 3^^ — 1) arctanf-u + x) — arctanf-u — x) 
X L 

+ -{x^ - 3u^x -3x- 2u^) ln[l + {u + xf] 



+ - 3m\ - 3x + 2u^) ln[l + (m - , 
with the coefficient r defined after eq.(ll). The associated three-body term reads 



(13) 



(70 + lr)9\Mmt [^^^ _ _ + flx^ _ 31^2 _ 5 _ (,2 _ ^2^2 



(87r)3/4 

3 ln(l -|- 4x^) -|- ( 6x —] arctan2x -|- arctan(-u + x) — arctan(-u — x 



2x 



2x 



{u' - x'Y + (12m^ + 27) (m^ - x^) + 
2u^^ 



+ + - 3x2 ^ ^ j 2^^^ + (m + x)2 
+ (^6 + - 3x2 _ ^ j j^^^ + (m - x)2j ^ , 
and the four-body term is given by the expression: 



(14) 



(8vr)V.^ 



-|-61n(l -|- 4x ) -|- ( 12x ) arctan2x -|- arctan(u + x) — arctan(u — x) 



2x 



y - u^f + 3x^ + Qu^x^ - 9u^ + 27x2 - ISu^ - 7 



+2(x2 - 3^2 - 3) In {[I + (m + x)2][l + (m - x)2]) + {u^ - x^f 



(15) 



The additional contributions to W{p^kf) from the iterated Ivr-exchange Hartree and Fock 
diagrams are collected in eqs. (20-25) of ref. [HI. In the appendix we present also novel analytical 
expressions for the three-body and four-body Hartree term of iterated Ivr-exchange. The total 
imaginary single-particle potential evaluated at zero nucleon momentum {p = 0) can even be 
written clS 8b closed form expression: 



WiO,kf) 



37^g\Mmi 
(4vr/.)4 



3 + 2(7o2 + 7? - 70 - 7i)J ^ + (70 + 71 - 3)^2 



2^2 7r2 
1 + m2 + 12 



+Li2(-1 -^2) + 



5 - 70 - 71 + ln(2 + u^)-- ln(l + u'^) 



ln(l + M2)l (16) 



where Li2(— a^^) = /J^ (iC (C + '3')~^lnC denotes the conventional dilogarithmic function. Note 
that eq.(16) determines the density dependence of the half-width of nucleon hole states at the 
bottom of the Fermi sea. 
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Figure 8: The imaginary single-particle potential W{p,kfQ) at saturation density versus the 
nucleon momentum p. The full curve corresponds to the parameter set constrained by the 
neutron matter equation of state. The dashed curve shows the case of maximal saturation 
density p^"""" = 0.191 fm'l 



Fig. 8 shows the imaginary single-particle potential W{p,kfo) at saturation density as a 
function of the nucleon momentum p. The full curve corresponds to the parameter set con- 
strained by the neutron matter equation of state (where fc/o = 248 MeV). The resulting value 
W{0, kfo) = 43.3 MeV is somewhat larger than the value W{0, kfo) ~ 40 MeV obtained in the 
self-consistent Brueckner calculation of ref . |2ZI using the phenomenological Paris NN-potential. 
The dashed curve in Fig. 8 shows the case of maximal saturation density p^^^ = 0.191 fm~^. 
With associated parameters 7o = 7i = 0.39 the imaginary single-particle potential at p = 
comes out as W{0, kfo) = 22.7 MeV. This value is close to the result W{0, kfo) ~ 20 MeV found 
in the calculation of ref. [23 employing the Gogny Dl effective interaction. For comparison the 
value W{0, kfo) = 29.7 MeV jH] results from the pure pion-exchange dynamics. The momentum 
dependence of W{p, kfo) is generic. As a consequence of the decreasing phase space available 
for redistributing the nucleon-hole state energy the curves in Fig. 8 drop with momentum p and 
W{p, kfo) reaches zero at the Fermi surface p = kjQ. According to Luttinger's theorem [SU^ this 
vanishing is of quadratic order, W{p,kf) ~ — kf)"^, a feature which is clearly exhibited by 
the curves in Fig. 8. 



6 Asymmetry energy 

This section is devoted to the density dependent asymmetry energy A{kf). The asymmetry 
energy is generally defined by the expansion of the energy per particle of isospin-asymmetric 
nuclear matter (described by different proton and neutron Fermi momenta kp^n = ^/(l -F <^)^''^) 
around the symmetry line: Eas{kp,kn) = E{kf) + 5'^ A{kf) + 0{5^). Evaluating the first dia- 
gram in Fig. 1 and combining that partial result with all cut-off dependent terms from iterated 
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diagrams one obtains the following contribution (linear in density) to the asymmetry energy: 



Aikf) 



3(477/, 



^(sro-Ti; 



:i7) 



with the parameters Fq and Fi defined in eqs.(2,5). From the second pion-exchange diagram 
in Fig. 1 with two medium insertions on parallel nucleon lines one gets a (finite) contribution 
to the asymmetry energy of the form: 



g\Mm't I 2 
— 2(7o + 7i)m + I670M arctan2M 



3(87r)3/4 



(670 + 27i)m + 



7o + 7i 
2u 



ln(l + 4m^ 



(18) 



The same (pion-exchange) diagram with three medium insertions gives rise to a contribution 
to the asjTiimetry energy which can be represented as a double-integral: 



A{kf) 



uxy{llu'' — ISx^y^^ 
3(^2 — x^y^) 



+ (u^ — 5x^?/^ 



H 



ln(l + 



+ 



X 



3(1 + S2) 6(1 + S2)2 

8s(l + s'){5s' - 3s - s") + (3 + s2)(8ss' - 35^ - 8s 



2uxy + {u — X y )H 



I 2^ 



+47iM^ 



2-ux?/[ln(l + 



3(^2 — x'^y'^) 



+ (^ln(l 



3(1 + 5^ 



(19) 



with s' = uds/du and s" = u'd'^s/du' denoting partial derivatives. In the case of the third 
diagram in Fig. 1 with three medium insertions the occurring double-integral can actually be 
solved. The /c^-contribution to the asymmetry energy originating from the iterated NN-contact 
interaction reads: 



7?(3-ln2)- 70^(2 + ln2) 



(20) 



The absence of an interference term proportional to 7071 is highly remarkable. The asymmetry 
energy A{kf) is completed by adding to the terms in eqs. (17-20) the contributions from the 
(relativistically improved) kinetic energy, (static) Ivr-exchange and iterated Ivr-exchange written 
down in eqs. (20-26) of ref.jH]. 

In Fig. 9, we show the asymmetry energy A{kf) versus the nucleon density p = 2/i;j/37r^. 
The full curve corresponds to the parameter set constrained by the neutron matter equation of 
state. The predicted value of the asymmetry energy at saturation density po = 0.134fm~'^ is 
A{kfo) = 30.8 MeV. This is close to the empirical value A(fc/o) — 33 MeV of ref.|[TH| obtained in 
extensive fits to nuclide masses. As a consequence of the built in constraints from the neutron 
matter equation of state the full curve in Fig. 9 starts to bend downwards only above relatively 
high densities p > 0.3 fm~^. The dashed curve in Fig. 9 shows the case of maximal saturation 
density p™^^ = 0.191 fm~^ realized by the parameters 7o = 7i = 0.39. We use furthermore 
the decomposition^ Fq + Fi = 6.04 + 0.42 as it follows from choosing a common cut-off Aq = 
Ai = 0.60 GeV. The resulting value A{kfQ) = 35.0 MeV of the asymmetry energy at saturation 

^Returning to eq.(2) one notices that Fi and 71 are necessarily approximately equal if 71 « 1/2. 
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Figure 9: The asymmetry energy A{kf) versus the nucleon density p = 2A;j/37r^. The full curve 
corresponds to the parameter set constrained by the neutron matter equation of state. The 
dashed curve shows the case of maximal saturation density p^^^ = 0.191 fm~^. 

density p^^^ = 0.191 fm~^ is also in good agreement with the empirical value. However, in 
that case the downward bending of A{kf) sets in immediately above saturation density. The 
situation concerning the high density behavior of the asymmetry energy is still controversial. 
While Brueckner-Hartree-Fock calculations of ref.jSI] predict a monotonic increase with density 
the results of the Urbana group show a downward bending of A{kf) at densities p = 
0.4 fm~^ or higher (see Fig. 8 therein). Altogether, it seems that the value of the asymmetry at 
saturation density is rather stable while the detailed density dependence of A{kf) reflects more 
pronouncedly some properties of the corresponding neutron matter equation of state. 



7 Nuclear matter at finite temperatures 

In this section, we discuss nuclear matter at finite temperatures T < 30 MeV. We are particu- 
larly interested in the first-order liquid-gas phase transition of isospin-symmetric nuclear matter 
and its associated critical point {pc,Tc). As outlined in ref.[13j a thermodynamically consistent 
extension of the present (perturbative) calculational scheme to finite temperatures is to relate 
it directly to the free energy per particle F{p, T), whose natural thermodynamical variables are 
the nucleon density p and the temperature T. In that case the free energy density pF{p, T) 
of isospin-symmetric nuclear matter consists of a sum of convolution integrals over interaction 
kernels JCj multiplied by powers of the density of nucleon states in momentum space: 



dip,) - 



27r2 



p]-2Mp 
1 + exp — — — — — 



(21) 



The effective one-body "chemical potential" p{p,T) entering the Fermi-Dirac distribution in 
eq.(21) is determined by the relation to the particle density p = 4 dpi pid{pi). We now write 
down the additional interaction kernels arising from the NN-contact interaction proportional 
to 7o,i. The contact interaction at first order and the cut-off dependent pieces from iterated 
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diagrams are summed up via the parameters Fq,! 
reads: 



Their contribution to the two-body kernel 



/C 



(r) 



{To + Ti)piP2 



(22) 



Obviously, this kernel generates a temperature independent contribution to the free energy 
per particle which just grows linear in density, F{p,T) ~ p. From the second pion-exchange 
diagram in Fig. 1 one derives a further (finite) contribution to the two-body kernel of the form: 



647r/4 



In 



ml + (pi + 
ml + (pi - P2)' 



-2{pi + P2) arctan "^^^ + 2{pi 



P2) arctan 



Pi -P2 



m-n 



(23) 



Temperature and density Pauli blocking effects in intermediate NN-states are incorporated in 
the three-body kernel /C3. The second and third diagram in Fig. 1 together give rise to the 
following contribution to the three-body kernel: 



/C 



(7) 



P1+P2 



dq 



|Pl~P2 



2 I 2 

7o +7i 



(70 + 7i)g^ 
ml + g2 



In 



\pI -pl + q^ + 2p3g| 



\pI 



pl + q^ 



2p3g| ■ 



(24) 



The remaining kernels building up the free nucleon gas part and the interaction contributions 
from Ivr-exchange and iterated Ivr-exchange are written down in eqs. (4,5,6,7,11, 12) of ref.|13j. 
Again, the (higher order) relativistic 1/M^-correction to the Ivr-exchange two-body kernel k!^2^^ 
(see eq.(5) in ref. JH]) is neglected here for reasons of consistency. A special feature of the finite 
temperature formalism is the so-called anomalous contribution pA{p, T) with no counterpart in 
the calculation to the ground state energy density pE{kf) at T = j32||33j. The basic physical 
mechanism behind this "anomalous" contribution is the possibility of interactions between 
particles and holes of the same momentum at finite temperature. From the combination of pion- 
exchange and zero-range NN-contact interaction we derive the following anomalous contribution 
to the free energy per particle of isospin-symmetric nuclear matter: 



Mp. t) 



2pfi'o'(p,T) 

A 



(70 + 71 



dpi I dp2 I dp3 d{pi)d{p2)[27r d{p2) - P2]d{p3) 



m: 



ml + {pi + P2) 



l)p,+ --^\n 

Ap2 ml + [pi 



P2)\ 

X 70 + 71 - l P3 + ^ In ^ , 

4p2 rnl + {p3-p2)^ 



(25) 



with the /i-derivative of the grand canonical potential due to static Ivr-exchange and the NN- 
contact couphng: 



n[{p,T) 



(pi 



dpi 

P2f 



dp2 d{pi 



d{p2 



{P1+P2) 



ml + [pi -P2Y 
and the second /i-derivative of the free nucleon gas part: 

poo 

nl{p,T) = -AM / dp 



P2 [ml + {pi + P2) 
-2(7o + 7i)pi 



d{pi) 
Pi 



(26) 



(27) 
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Figure 10: Pressure isotherms P{p,T) of isospin-symmetric nuclear matter for the case of 
maximal saturation density. The coexistence region of the liquid and gas phase ends at the 
critical point pc — 0.09 fm~^ and Tc ~ 21 MeV. 

The first term in eq.(25) originates from taking into account the contributions of (static) Ivr- 
exchange and the NN-contact coupling in the Legendre transformation from the grand canonical 
potential to the free energy density and from the perturbative shift of the chemical potential 
/i ^ /i — Q[{p,T) /Qq{p,T) (for details on that procedure, see ref.p3|). We have explicitly 
checked that the anomalous contribution vanishes identically at T = for all densities p. This 
vanishing, A{p, 0) = 0, is an automatic consequence of the Kohn-Luttinger-Ward theorem 
j^n] which of course hold in our case since the Fermi surface and the NN-interactions are 
invariant under spatial rotations. Furthermore, it is interesting to observe that the temperature 
and density dependent anomalous contribution A{p, T) vanishes identically in the chiral limit 
m^r = 0. The latter feature has the practical consequence that the behavior of nuclear matter 
at finite temperatures T < 30 MeV is only marginally influenced by the anomalous contribution 
A{p,T). 

Fig. 10 shows for the case of maximal saturation density (70 = 71 = 0.39) the calculated 
pressure isotherms P{p,T) = p'^dF[p,T)/dp of isospin-symmetric nuclear matter at four se- 
lected temperatures T = 0, 15, 21, 25 MeV. As it should be these curves display a first-order 
liquid-gas phase transition similar to that of the van-der-Waals gas. The coexistence region 
between the liquid and the gas phase (which has to be determined by the Maxwell construction) 
terminates at the critical temperature T^. From there on the pressure isotherms P{p,T) grow 
monotonically with the nucleon density p. We find here a critical temperature of Tc ~ 21 MeV. 
This is already compatible with the value Tc = (20 ± 3) MeV as extracted in ref. [S3 from 
multi-fragmentation data in proton on gold collisions. Most other calculations get still lower 
values of the critical temperature, typically Tc ~ 18 MeV [13 ESI EE]- The reduction of Tc in 
comparison to Tc ~ 25.5 MeV obtained in ref. [IS] from pion-exchange dynamics alone results 
from the less strong momentum dependence of the real single-particle potential U{p, kfo) near 
the Fermi surface p = fc/o (see Fig. 7). As a general rule Tc grows with the effective nucleon mass 
M*{kfo) at the Fermi surface. Because of that property Tc ~ 21 MeV is the absolute minimum 
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of the critical temperature possible within our perturbative calculation to fourth order. 

For the parameter set (70 = 0.39, 71 = —0.75) constrained by the equation of state of pure 
neutron matter the unrealistic momentum dependence of the associated real single-particle 
potential U{p, kfo) (see Fig. 6) ruins the behavior of nuclear matter at finite temperatures. The 
critical temperature of the liquid-gas phase transitions would lie in that case around Tc ~ 
40 MeV, thus exceeding realistic values by more than a factor of two. 

Evidently, the single-particle properties around the Fermi surface are decisive for the ther- 
mal excitations and therefore they crucially influence the low temperature behavior of nuclear 
matter. 



8 Nuclear energy density functional 



The energy density functional is a general starting point for nuclear structure calculations 
within the framework of the self-consistent mean-field approximation ^7\. In this context 
effective Skyrme forces jHHl (HHl HO] have gained much popularity because of their analytical 
simplicity and their ability to reproduce nuclear properties over the whole periodic table. In a 
recent work |41j, we have calculated the nuclear energy density functional which emerges from 
(leading and next-to-leading order) chiral pion-nucleon dynamics. The calculation in ref. [H] 
included (only) the Ivr-exchange Fock diagram and the iterated Ivr-exchange Hartree and Fock 
diagrams. Therefore, the interest here is on the additional contributions from the NN-contact 
interaction appearing at the same order in the small momentum expansion. 

Going up to quadratic order in spatial gradients (i.e. deviations from homogeneity) the 
energy density functional relevant for N = Z even-even nuclei reads 



£[p,T]=pEikf) + 



T 



-Pk) 



1 5fc^ 
2M 56M3 



+ Fr{kf) 



+ {VpfFy{kf) 



(28) 



with p(r ) = 2A;j(r )/37r^ the local nucleon density (expressed in terms of a local Fermi momen- 
tum kf{f)) and T{f) the local kinetic energy density. Note that we have left out in eq.(28) the 
spin-orbit coupling term since the NN-contact vertex proportional to 70.1 does not generate any 
contribution to it. The starting point for the construction of an explicit nuclear energy density 
functional S[p,t] is the bilocal density matrix as given by a sum over the occupied energy 
eigenfunctions: J2aeocc'^ai^ ~ fl/2)\E'j^(r + a/2). According to Negele and Vautherin jl2j it 
can be expanded in relative and center-of-mass coordinate, a and f, with expansion coefficients 
determined by purely local quantities (nucleon density, kinetic energy density and spin-orbit 
density). As outlined in sect. 2 of ref. [H] the Fourier transform of the (so-expanded) density 
matrix defines in momentum space a medium insertion for the inhomogeneous many-nucleon 
system: 



d^r e-"'-' e{kf - \p\)\ 1 + 



357r2 
8k} 



(5p' 



3k') 



(29) 



In comparison to eq.(5) in ref. |^ we have left out here the spin-dependent part (generating 
spin-orbit couplings) since it does not come into play. The strength function F^{kf) in eq.(28) 
emerges via a perturbation on top of the density of states 6{kf — |p |). As a consequence of that 
Frikf) can be directly expressed in terms of the real single-particle potential U{p^ kf) as: 



Fr{kf) 



35 
4^ 



dpp'i^p' — 3A;?) U{p, kf 



(30) 
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Note that any p-independent contribution, in particular, the cut-off dependent ro,i-term in 
eq.(9) drops out. In the medium insertion eq.(29) r — Spkj/d is accompanied by — V^p/4. 
After performing a partial integration one is lead to the following decomposition: 



Mkf) 



8k} 



dki 



+ Ukf) 



(31) 



where Fct{kf) comprises all those contributions for which the (Vp)^-factor originates directly 
from the momentum-dependence of the interactions. We enumerate now the contributions to 
the strength functions Fr^d{kf) generated by the NN-contact vertex. From the second pion- 
exchange diagram in Fig. 1 with two medium insertions we get: 
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while the same diagram with three medium insertions gives: 
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In case of the last diagram in Fig. 1 with three medium insertions we can even solve the integrals 
and get: 

Evidently, there is no contribution to Fct{kf) from this last diagram since the NN-contact 
vertex is momentum independent. The additional contributions to F^^dikf) from Ivr-exchange 
and iterated Ivr-exchange can be found in sect. 3 of ref. 



The expression in eq.(28) multiplying the kinetic energy density T(r ) has the interpretation 
of a reciprocal density dependent effective nucleon mass: 

M*{p) = M 

In Fig. 11, we show the ratio effective over free nucleon mass M*{p)/M as a function of the 
nucleon density p = 2/c^/37r^. The dashed curve shows the case of maximal saturation density 



bk ~ 

1 ^— + 2MFJkf) 

28M2 ^' 



(37) 



19 



1.15 



1.05 




0.95 



0.85 



p[fm1 



Figure 11: The effective nucleon mass M*{p) divided by tlie free nucleon mass M versus tlie 
nucleon density p = 2fc^/37r^. Tlie daslied curve sliows tfie case of maximal saturation density 



Po 



0.191 fm"'^. The dotted curve gives the result of the pion-exchange contributions alone. 



^max _ Q_x91fm ^ while the dotted curve gives the result of the pion-exchange contributions 
alone jH]. As a consequence of the more moderate momentum dependence of the real single- 
particle potential U{p, kj) (see Fig. 7) the effective nucleon mass M*{p) is now reduced for all 
densities p < 3po/2. In comparison to most (non-relativistic) mean- field calculations [H ^ r ^ HH] 
the reduction by at most 15% as obtained here is a rather weak one. Note that the deviation of 
the effective nucleon mass M*{p) from the free nucleon mass M does not just grow linear with 
density, but there are very strong curvature effects. This has mainly to do with the explicit 
presence of the small mass scale m.^ = 135 MeV in our calculation. For the parameter set 
(70 = 0.39, 7i = —0.75) constrained by the neutron matter equation of state the unrealistic 
momentum dependence of the corresponding real single-particle potential U{p, kfo) (see Fig. 6) 
would even lead to a singularity of M*{p) at p ~ 0.16 fm~^. 



In Fig. 12, we show the strength function Fy{kf) related to the (Vp)^-term in the nuclear 
energy density functional versus the nucleon density p = 2fcj/37r^. The full curve corresponds 
to the parameter set constrained by the neutron matter equation of state. The dashed curve 
shows the case of maximal saturation density p^^^ = 0.191 fm~^. For comparison, the dotted 
curve gives the result of the pion-exchange contributions alone jH]. The three horizontal lines 
represent the constant values F\r{kf) = [9ti — (5 + 4a;2)i2]/64 of the Skyrme forces Sly pi] . 
Sill and MSk [lOj. In the case of Sly and MSk we have performed averages over the 
various parameter sets Sly 4-7 [22] and MSkl-6 [40J. One observes that in the case of maximal 
saturation density (dashed line in Fig. 12) the effects from the NN-contact interaction move the 
prediction for F^{kf) further away from the empirical values. For the parameter set constrained 
by the neutron matter equation of state one gets from the NN-contact interaction a substantial 
lowering of the strength function Fsj{kf) compared to the pure pion-exchange contributions 
[llj . In the region pq/2 < p < 3po/2 the prediction falls now into the empirically allowed band. 
The strong rise of Fv(/i;/) at very low densities has to do with the presence of the small mass 
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Figure 12: The strength function Fy{kf) related to the (Vp)^-term in the nuclear energy 
density functional versus the nucleon density p = 2kj^/3n'^. The full curve corresponds to the 
parameter set constrained by the neutron matter equation of state. The dashed curve shows 
the case of maximal saturation density p™^ = 0.191 fm~^. The dotted curve gives the result 
of the pion-exchange contributions alone. 

scale = 135 MeV, and with associated chiral singularities (of the form and m~^). 

Altogether, one finds a conflicting situation concerning the effects of the NN-contact inter- 
action on the strength functions Fr{kf) and F^{kf). The parameter set with an unrealistic 
p-dependencc of single-particle potential U{p, kf), and therefore an unrealistic effective nucleon 
mass M*{p), leads to an improved description of the strength function Fy{kf ) related to the 
(Vp)^-term in the nuclear energy density functional. The inclusion of short-range NN-terms 
does therefore not lead to a consistent improvement of inhomogeneous nuclear matter proper- 
ties. 

9 Summary and conclusions 

In this work we have extended a recent chiral approach to nuclear matter by including the 
effects from the most general (momentum- independent) NN-contact interaction. By iterating 
the combination of one-pion exchange and this two-parameter contact-vertex to second order 
the emerging energy per particle exhausts all terms possible up-to-and-including fourth order in 
the small momentum expansion. Two cut-offs Ao,i, depending on the total NN-isospin / = 0, 1, 
are introduced to regularize the (primitive) linear divergences dl 1 of some three-loop in- 
medium diagrams. Effectively, the cut-offs Ao,i parameterize the strength of a term linear in 
density, see eqs.(2,5). The equation of state of pure neutron matter, En{kn), can be reproduced 
very well up to quite high neutron densities of p„ = 0.5 fm~^ by adjusting the strength 71 of a 
repulsive nn-contact interaction. With 71 = —0.75 it is comparably strong as the Ivr-exchange. 

Binding and saturation of isospin-symmetric nuclear matter is a generic feature of the 
present perturbative calculation. Fixing the maximum binding energy per particle to —E{kfo) 
— 15.3 MeV we find that any possible equihbrium density po lies below p^'^"^ — 0.191 fm~^. The 
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constraint 71 = —0.75 from the neutron matter equation of state leads, however, to a somewhat 
too low saturation density of po = 0.134 fm^^. Fixing on the other hand the saturation density 
to Po = 0.158 fm"^ 1221 we find that the maximum binding energy per particle —E{kfQ) lies 
above 10.7 MeV. 

We have also investigated the effects of the NN-contact interaction on the complex single- 
particle potential U{p, kf) + i W{p, kf). We have found that the effective nucleon mass at the 
Fermi-surface is bounded from below by M*{kfo) > 1AM. The minimal effective nucleon 
mass occurs simultaneously together with the maximal saturation density. The corresponding 
parameter set 70 = 7i = 0.39 describes the situation in which the explicitly introduced short- 
range NN-terms cancel almost exactly the contact pieces of l7r-exchange. The parameter set 
constrained by the neutron matter equation of state leads to unrealistic real and imaginary 
single-particle potentials, U{p, kfo) and W{p, fc/o), as in the scheme of Lutz et al. [3 Ell. The 
downward bending of the asymmetry energy A{kf) above nuclear matter saturation density is 
a generic feature of our perturbative (fourth-order) calculation. The value of the asymmetry 
energy at saturation density A{kfQ) is rather stable and comes out close to its empirical value 
A{kfo) ~ 33 MeV. 

The properties of nuclear matter at finite temperatures have also been investigated. The 
critical temperature of the liquid-gas phase transition has a minimal value of Tc ~ 21 MeV. 
Good single-particle properties around the Fermi surface are a necessary condition for a realistic 
behavior of isospin-symmetric nuclear matter at finite temperatures T < 40 MeV^ 

The influence of the NN-contact interaction on the effective nucleon mass M*(p) and the 
strength function Fy{kf) of the (Vp)^-term in the nuclear energy density functional S[p,t] 
have also been studied. For the parameter set with reduced M*{p) the prediction for Fy{kf) 
is moved further away from the empirical band. On the other hand the parameter set with 
unrealistic M*(p) substantially improves F\/{kf). 

Altogether, there is within this complete fourth order calculation no "magic" set of short- 
range parameters with which one could reproduce simultaneously and accurately all semi- 
empirical properties of nuclear matter. The conditions for a good neutron matter equation of 
state and for good single-particle properties (and consequently a realistic finite temperature 
behavior) are in fact mutually exclusive. This clearly points at inherent limitations of the 
perturbative chiral approach truncated at fourth order in the small momentum expansion. 
While complete calculation at fifth order or even sixth order are rather demanding (four-loop 
and five-loop diagrams need to be evaluated) one might also combine the chiral approach (to 
nuclear matter) with non-perturbative methods as it is done commonly for free NN-scattering 

1201 EI!- 

Appendix 

In this appendix we take the opportunity to write down closed form analytical expressions 
for the three-body and four body contributions W^^^lp, kf) to the imaginary single-particle 
potential as they arise from the iterated Ivr-exchange Hartree diagram. Solving the integrals in 
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eq.(22) of ref-jH] gives for the three-body term: 



/21 
2x 

1 



35 + IhW + 25m^ - 



197a;2 



- 2u'x' - 7x^ + 18 ln(l + Ax^) 



+ 1 30a; 1 arctan 2x + 



arctanfu + x) — arctanfu — x) 



2x 



{x^ - u^f + 9(4^2 + 15) (x^ - M^) - 56 
2u^ 

X 



+a(3x'^ 9«^) ln[l + {u + xf] 

+4 (3x2 _ 9 ^ ^ _ j ^ _ ^)2] I 



(38) 



with X = p/m.^ and u = kf/mT^. After complete evaluation of its integral representation given 
in eq.(24) of ref.pj the four-body Hartree term reads: 



W,{p,kf) 



"J -13u^- 7- 2^2(49 + 3x2) + — 361n(l + 4x2) 



arctan(M + x) — arctan(n — x) 



/ 21 

+ 1 60x 1 arctan 2x + 

\ X 

1 



X ■ 



2x 



- x^f - 9x^ - ISu^x^ + 27u^ - 135x2 + 75^^ + 49 



+4(9 + - 2x2) In {[I + {u + xf] [I + {u - xf]^ | . 



(39) 



Combining these novel results with the two-body Hartree term W2{p, kj) written in eq.(20) of 
ref.jnj one immediately verifies Luttinger theorem PU] : 



W^ikf, kf) + W^{kf, kf) + W,{kf, kf) = 0, 



and the zero sum-rule jH]: 



JO 



dpp^ 



\W2{P, kf) + ^Wsip, kf) + \w,{p, kf) 



0. 



(40) 



(41) 



The latter statement expresses the fact that the energy per particle E{kf) can be reconstructed 
from the two-, three- and four-body contributions to the single-particle potential (graphically 
this is done by closing the open nucleon line of self-energy diagrams) and that no imaginary 
part is left over from this procedure. 
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